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9. = 1 — 2P has a stable equilibrium at P = = —2(1 — 2P) 




10. 2P) has an unstable equilibrium at P — 0 and a stable equilibrium at P = 1. 

§ = (l-4P)f =P(1-4P)(1-2P) 
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11. Yjy — 2P(P - 3) has a stable equilibrium tit P — 0 and an unstable equilibrium at P = 3. 
§ = 2(2P - 3)f = 4P(2P - 3)(P - 3) 



p '> o ! p'<o ! f>o 
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12. = 3P( 1 — P) (P — |) has a stable equilibria at P = 0 and P — 1 an unstable equilibrium at P = 

W = -K 6p2 - 6P+1 )f = l p ( p - ( p - 1) ( p - 3 ~¥ 1 ) ( p - 1) 
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13. 

Before Catastrophe After Catastrophe 



P P 




Before the catastrophe, the population exhibits logistic growth and P(t) — > Mo, the stable equilibrium. After the 
catastrophe, the population declines logistically and P(t) — > Mi, the new stable equilibrium. 

14. § = rP(M - P)(P - m), r, M, m > 0 
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The model has 3 equilibrium points. The rest point P = 0, P = M are asymptotically stable while P = m is unstable. For 
initial populations greater than m, the model predicts P approaches M for large t. For initial populations less than m, the 
model predicts extinction. Points of inflection occur at P = a and P = b where a — I[m) m — \J M 2 — mM + m 2 ] and 
b = 5 [m + m + \/M 2 — mM + m 2 ] • 

(a) The model is reasonable in the sense that if P < m, then P — > 0 as t — > oo; if m < P < M, then P — > M as t - - oo; if 
P > M, then P — > M as t — > oo. 

(b) It is different if the population falls below m, for then P — > 0 as t — > oo (extinction). If is probably a more realistic 
model for that reason because we know some populations have become extinct after the population level became too 
low. 
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(c) For P > M we see that ^ = rP(M P)(P m) is negative. Thus the curve is everywhere decreasing. Moreover, 

P = M is a solution to the differential equation. Since the equation satisfies the existence and uniqueness conditions, 
solution trajectories cannot cross. Thus, P — > M as t — > oo. 

(d) See the initial discussion above. 

(e) See the initial discussion above. 



15. ^ = g — jA 2 , g, k, m > 0 and v(t) > 0 
Equilibrium: ^ = g 



k v 2 = 0=>v= ,/f 



Concavity: g = -2(£v) | = -2(£v) (g - ^v 2 ) 



(a) 



(b) 



o A, 

dt 2 ' 



“• > i < d ' 



. dt 2 ■ 



mg 

u eq=V k 




(C) Vtermina! = \/ rfm = 178.9 f = 122 mph 



16. 



F = Fp - F r 
ma = mg — 

S=g-^.v(0)=v 0 
Thus, ^ = 0 implies v = (™ g 
terminal velocity; if Vq > 



)”, the terminal velocity. If Vo <m , the object will fall faster and faster, approaching the 
)", the object will slow down to the terminal velocity. 



17. F = F p - F r 
ma = 50 — 5|v| 

S = i(50-5|v|) 

The maximum velocity occurs when ^ = 0orv = 1 0 ^ . 



18. (a) The model seems reasonable because the rate of spread of a piece of information, an innovation, or a cultural fad is 
proportional to the product of the number of individuals who have it (X) and those who do not (N — X). When X is 
small, there are only a few individuals to spread the item so the rate of spread is slow. On the other hand, when 
(N — X) is small the rate of spread will be slow because there are only a few indiciduals who can receive it during the 
interval of time. The rate of spread will be fastest when both X and (N — X) are large because then there are a lot of 
individuals to spread the item and a lot of individuals to receive it. 

(b) There is a stable equilibrium at X = N and an unstable equilibrium at X = 0. 

^ — k^jjr (N — X) — kX^ = k 2 X(N - X)(N - 2X) => inflection points at X = 0, X = f , and X = N. 
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(d) The spread rate is most rapid when x = ^ . Eventually all of the people will receive the item. 

19. L| + Ri = V^| = X_Ri = R(v_i),v jL ,R>0 
Equilibrium: | = |(|-i)=0=M=| 

Concavity: § = -(!)! = ~ 0 

Phase Line: 




. _ v_ 

'eq R 

If the switch is closed at t = 0, then i(0) = 0, and the graph of the solution looks like this: 




As t — > oo, it — i ► i s teady state = p- (In the steady state condition, the self-inductance acts like a simple wire connector and, as 
a result, the current throught the resistor can be calculated using the familiar version of Ohm's Law.) 

20. (a) Free body diagram of the pearl: 




W -mg 



(b) Use Newton's Second Law, summing forces in the direction of the acceleration: 
mg - Pg - kv = ma => f = ( m 5r f )g - £v. 
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(d) 

V 




9.5 SYSTEMS OF EQUATIONS AND PHASE PLANES 

1. Seasonal variations, nonconformity of the environments, effects of other interactions, unexpected disasters, etc. 

2. x = r cos # => = -rsin# 1 ^ + cos#| = y + x - x(x 2 + y 2 ) = rsin# + rcos# - Ecos# 

y = rsin# => ^ = rcos# ^ + sin# g 2 = — x + y — x(x 2 + y 2 ) = — rcos# + rsin# — Esin# 

Solve for ^ by adding cos # x eq(l) to sin# x eq(2): 

cos 2 6 gj + sin 2 9 $ = cos#(rsin# + rcos# — Ecos#) + sin#(— rcos# + rsin# — Esin#) 

=> — rsin#cos# + rcos 2 9 — r 3 cos 2 9 — rsin#cos# + rsin# — r 3 sin 2 9 = r — r 3 = r(l — r 2 ) 

Solve for ^ by adding (— sin#) x eq(l) to cos 9 x eq(2): 

rsin 2 #^ + rcos 2 9 ^ = — sin#(rsin# + rcos# — Ecos#) + cos#(— rcos# + rsin# — r 3 sin#) 

=> r^ = —rsin 2 # — rsin# cos# + Esin# cos# — rcos 2 # + rsin# cos# — Esin# cos# = — r =£■ ^ = — 1 

= (1)^1 — (l) 2 j = 0, thus the trajectory 

remains on the circle, and rotates around the circle in a clockwise direction, since ^ = — 1. The solution is periodic since 
at any point (x, y) on the trajectory, (x, y ) = (r cos #, r sin #) = (1 cos #, 1 sin #) = (cos #, sin #) => both x and y are 
periodic. 

3. This model assumes that the number of interactions is porportional to the product of x and y: 

| = (a-by)x, a < 0, f = m(l - ^)y-nxy = y(m- f 2 y -nx). 

To find the equilibrium points: 

(h = 0=^ (a - b y)x = 0 => x = 0 or y = jj 

(remember g < 0); ^ °’ M ^ 

!= 0 ^y(m-Sy_ nx ) ^ y = Oory = -^fx + M; 

Thus there are two equlibrium points, both occur when 
x = 0, (0, 0) and (0,M). 



(o. 




If r = 1 (that is, the trajectory starts on the circle x 2 + y 2 — 1 ), then ^ 
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Implies coexistence is not possible because eventually trout 
die out and bass reach their population limit. 




4. The coefficients a, b, m, and n need to be determined by sampling or by analyzing historical data. Then, more specific 
graphical predictions can be made. These predictions would then have to be compared to actual population growth 
patterns. If the predictions match actual results, we have partially validated our model. If necessary, more tests could be 
run. However, it should be remembered that the primary purpose of a graphical analysis is to analyze the behavior 
qualitatively. With reference to Figure 9.29, attempt to maintain the fish populations in Region B through stocking and 
regulation (open and closed seasons). For example, should Regions A or D be entered, restocking the appropriate species 
can cause a return to Region B. 

5. (a) Logistic growth occurs in the absence of the competitor, and simple interaction of the species: growth dominates the 

competition when either population is small so it is difficult to drive either species to extinction. 

(b) a = per capita growth rate for trout 
m = per capita growth rate for bass 

b = intensity of competition to the trout 
n = intensity of competition to the bass 
ki = environmental carrying capacity for the trout 
k 2 = environmental carrying capacity for the bass 

(c) ^ = 0 => a^l — j^x — bxy = a^l-^-j-byx=04>x = 0ora^l-|-j-by = 0=^x = 0or 

y=§-X x; *= 0=>m ( 1 ^£) y ~ nxy = [ m (' ~ i) - nx ]y = o=> y = oor 

m^l — — nx = 0=>y = 0ory = k 2 — There are five cases to consider. 

Case I: £ > k 2 and ^ > ki . 

By picking g > k 2 and ” > ki we ensure an equilibrium point exists inside the first quadrant. 




Graphical analysis implies four equilibrium points exist: (0, 0), (ki, 0), (0, k 2 ), and TmXnkX ) 

(the point of intersection of the two boundaries in the first quadrant). All of these equilibrium points are unstable 
except for the point of intersection. The possibility of coexistence is predicted by this model. 
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Case II: g > k 2 and f < ki . 
(0, k 2 ): unstable 
(ki,0): stable 
(0,0): unstable 
Trout wins: (ki, 0) 

Not sensitive 
No coexistence 



Case III: “ < k 2 and ® > kj . 
(0, k 2 ): stable 
(ki,0): unstable 
(0,0): unstable 
Bass wins: (0, k 2 ) 

Not sensitive 
No coexistence 



Case IV: £ < k 2 and 31 < ki . 

(0, k 2 ): stable 
(ki,0): stable 
(0,0): unstable 

( amkl ~ h b TV k C amk2 ~K a,1 , k ( k2 V unstable 

\ am — bnkik2 am — bnkjk2 ) 

Bass or trout: (0, k 2 ) or (ki, 0) 

Very sensitive 

Coexistence is possible but not predicted 





If we assume g < k 2 and “ < kj then graphical analysis implies four equilibrium poins exist: (0, k 2 ), (ki, 0), 

(0, 0), and ^ a am-^nkik^ 2 1 ^m-'bnkfe 2 ) 9he P°i nt °f intersection of the two boundaries in the first quadrant). 
Case V: g = k 2 and ^ (lines coincide). 

(0, k 2 ): stable 
(ki,0): stable 
(0,0): unstable 

Line segment joining (0, k 2 ) and (ki, 0): stable 
Bass wins: (0, k 2 ) 

Not sensitive 

Coexistence is likely outcome 



Note that all points on the line segment joining (0, k 2 ) and (ki, 0) are rest points. 




6. For a fixed price, as Q increases, gg gest smaller and, possibly, becomes negative. This observation implies that as the 
quantity supplied increases, the price will not rise as fast. If Q gets high enough, then the price will decrease. Next, 
consider For a fixed quantity, as P increases, gjggets larger. Thus, as the market price increases, the quantity 
supplied will increase at a faster rate. If P is too small, g2 will be negative and the quantity supplied will decrease. 
This observation is the traditional explanation of the effect of market price levels on the quantity supplied. 
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Now ^ > 0 when PQ < 20,000 and P > 0; ^ < 0 otherwise. ^ > 0 when P > ^ and Q > 0; ^ < 0 otherwise, 
(b) These considerations give the following graphical analysis: 




The equilibrium point (0, 0) is unstable. The graphical analysis for the point (25.8, 775) is inconclusive: trajectories 
near the point may be periodic, or may spiral toward or away from the point. 

(c) The curve ^ = 0 or PQ = 20000 can be thought of as the demand curve; ^ = 0 or Q = 30P can be viewed as the 
supply curve. 



7. (a) 

(b) 

(c) 



| = ax^bxy = (a-by)xand^ = my-nxy = (m-nx)y^ f = || =► % = f = 17=^ 

1 = 7^=* (f - b ) d y = (^ - n ) dx => /(^ - b )dy = /(f - n )dx aln |y| - hy = mln | X | nx + C 
=>■ ln|y a | + lne~ by = In | x m | + lne~ nx + lne c => ln|y a e~ by | = ln|x m e” nx e c | => y a e~ by = x m e“ nx e c , let K = e c 
=> y a e _by = Kx m e~ nx 

f(y) = y a e~ by => f '(y) = ay a ~ 1 e” by — by a e _by = y a ~ 1 e~ by (a — by) and f '(y) = 0=>y = 0ory = a ; 
f"(|) = — b(|) a *e _a < 0 => f(y) has a unique max of M y = (-^) d wheny=jj. g(x) = X m e _nx 
=> g'(x) = mx m_1 e” nx — nx m e~ nx = x m-1 e~ nx (m — nx) and g'(x) = 0=>x = 0orx = 
g"(™) = — n (^) m 'e - ™ < 0 => g(x) has a unique max of M x = (^) m whenx=™. 



f(y) g(x) 




(d) Consider trajectory (x, y) — > (”, |). y a e by = Kx m e nx => ^ ^ = K, taking the limit of both sides 
=> lim ( 4y ' ) = lim K => ^ = K. Thus, ^ represents the equation any solution 

x— >m/n \ e / x-»m/n x x 

y— >a/b y— >a/b 

trajectory must satisfy if the trajectory approaches the rest point asymptotically. 
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(e) Pick initial condition yo < 5 • Then, from the figure at 

Myx™ _ _y8 



= ^ < M y and thus 



right, f(yo) < M y implies 
^ < M x . From the figure for g(x), there exists a 
unique Xo < ™ satisfying x n ™ < M x . That is, for each 
y < g there is a unique x satisfying ^ ^ Thus, 

there can exist only one trajectory solution approaching 
(™, g) . (You can think of the point (xo, yo) as the initial 
condition for that trajectory.) 



y Bass 



-(•}- 



/ 

Unique Xq 



Trout 



(f) Likewise there exists a unique trajectory when yo > g- Again, f(yo) < M y implies ^ < My and thus 

^ < M x . From the figure for g(x), there exists a unique Xo > ™ satisfying x n '" < M x . That is, for eachy > g there is 
a unique x satisfying ^ = jy^ ,£• Thus, there can exist only one trajectory solution approaching (™, g). 



Le tz = y' = I 



dz 

dx 



= z' = y", then given the differential equation y" — F(x, y, y'), we can write it as the following 

A., 



system of first order differential equations: ^ = 



t = F(x, y, z) 



dx v ’ - - 

In general, for the n th order differential equation given by y (n ' = F(x, y, y', y", . . ., y (n_1) ), let z \ = y' = % 

fr = z \ = y". let z - = z i = y"» => g| = z 2 = y"'. • • - let z n-i = K~2 = y (n_ =>■ z^_j = y (n >. This gives us the 



ax 1 s 7 ^ 1 •/ 7 ax ^ 

following system of first order differential equations: ^ = Z] 

dzi 

dx 



= z 2 



^= Z 3 



= F(x, y, zi, z 2 , . . ., z n _ ! ) 

9. In the absence of foxes => b = 0 =£- gg = ax and the population of rabbits grows at a rate proportional to the number of 
rabbits. 

10. In the absence of rabbits =>d = 0=>^| = — cy and the population of foxes decays (since the foxes have no food source) 
at a rate proportional to the number of foxes. 

11 . gg = (a — by)x = 0 => y = g or x = 0; gg = (— c + dx)y = 0 => x = g ory = 0 => equilibrium points at ( 0 , 0 ) or 

(g, g) . For the point (0, 0), there are no rabbits and no foxes. It is an unstable equilibrium point, if there are no foxes, but 
a few rabbits are introduced, then gg = a => the rabbit population will grow exponentially away from ( 0 , 0 ) 

12. Let x(t) and y(t) both be positive and suppose that they satisfy the differential equations gg = (a — by)x and 

gg = (— c + dx)y. Let C(t) = alny(t) - by(t) - dx(t) + clnx(t) => C'(t) = a - by'(t) - dx'(t) + c ^ 

= (yW -b ) y, W + (l^j _d ) x '( t ) = (w) - b )(-c + dx(t))x(t)+ (^j — d) (a — by(t))y(t) = 0 
Since C'(t) = 0 => C(t) = constant. 

13. Consider a particular trajectory and suppose that (xo, yo) is such that xo < g and yo < g. then gg > 0 and g| < 0 =>- the 
rabbit population is increasing while the fox population is decreasing, points on the trajectory are moving down and to the 
right; if xo > g and yo < g, then gg > 0 and g? > 0 =>■ both the rabbit and fox populations are increasing, points on the 
trajectory are moving up and to the right; if xo > a and yo > g , then gg < 0 and ^ > 0 => the rabbit population is 
decreasing while the fox population is increasing, points on the trajectory are moving up and to the left; and finally if 

Copyright © 2010 Pearson Education, Inc. Publishing as Addison-Wesley. 




562 Chapter 9 First-Order Differential Equations 



xo < | and yo > g, then ^ < 0 and $ < 0 => both the rabbit and fox populations are decreasing, points on the trajectory 
are moving down and to the left. Thus, points travel around the trajectory in a counterclockwise direction. Note that we 
will follow the same trajectory if (xo, yo) starts at a different point on the trajectory. 

14. There are three possible cases: If the rabbit population begins (before the wolf) and ends (after the wolf) at a value larger 
than the equilibrium level of x = then the trajectory moves closer to the equilibrium and the maximum value of the 
foxes is smaller. If the rabbit population begins (before the wolf) and ends (after the wolf) at a value smaller than the 
equilibrium level of x = but greater than 0, then the trajectory moves further from the equilibrium and the maximum 
value of the foxes is greater. If the rabbit population begins and ends very near the equilibrium value, then the trajectory 
will stay near the equilibrium value, since it is a stable equilibrium, and the fox population will remain roughly the same. 



CHAPTER 9 PRACTICE EXERCISES 



1 . 



y' = xe y \/x 

=> -y = In 



— 2 => e y dy = x-y/x 

' -2(x-2f 2 (3x + 4) _ c j 



— 2 dx => — e 

=> y = -In [ 



— y 2(x — 2) 3/2 (3 x + 4) . p y 

15 -t- ^ e 

~2(x-2f 2 (3x + 4) _ c j 



— 2(x — 2) 3/2 (3 x 4- 4) r 

15 ^ 



2. y' = xye x2 => f = e x \dx => In y = ^ + C 

3. sec x dy + x cos 2 y dx = 0 => => tan y = —cos x — x sin x + C 

J J cos z y sec x J 

4. 2x 2 dx — 3y/ycsc x dy = 0 => 3y/ydy = ^^dx =>■ 2y 3/2 = 2(2 — x 2 )cos x + 4x sin x + C 
=> y 3/2 = (2 — x 2 )cos x + 2x sin x + Ci 

5- y' = | => ye~ y dy = f ^ (y + l)e“ y = -In |x| + C 

6. y' = xe x ~ y csc y => y' = ^csc y =>■ ( ^dy = xe x dx => ^(sin y — cos y) = (x — l)e x + C 

7. x(x - l)dy - y dx = 0 => x(x - l)dy = y dx =>• f = => In y = ln(x - 1) - ln(x) + C 

=S> In y = ln(x - 1) - ln(x) + In Ci => In y = ln( Cl(x x ~ 1} ) => y = 



y' = (y 2 - i)(x-‘) 



dy 

y 2 -i 



dx 



o=,„ x . 



ta (m) 



= 21n x + In Ci 



y- 1 
y+i 



= Cix 2 



9. 2y ' — y = xe x/2 => y' - \y = |e x/2 . 
p(x) = -i, v(x) = e/H) dx = e _x/2 . 
e _x/2 y' - |e- x/2 y = (e" x/2 ) (|) (e x/2 ) 

10. \ + y = e~ x sin x => y' + 2y = 2e _x sin x. 
p(x) = 2, v(x) = ef 2dx = e 2x . 

e 2x y' + 2e 2x y = 2e 2x e~ x sin x = 2e x sin x => gj(e 2x y) = 2e x sin x => e 2x y = e x (sin x — cos x) + C 
=>■ y = e~ x (sin x — cos x) + Ce~ 2x 



dx 



(e- x/2 y) = 



e x/2 y _ 



+ C^y = e x/2 (f +C) 
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11. xy' + 2y=l-x- 1 ^y'+(?)y=^^. 
v(x) = Q 2 If = e 21nx = e lnx2 = x 2 . 

x 2 y' + 2xy = x - 1 =>• |^(x 2 y) = x- l=>x 2 y=f-x + C=>y=l-? + g 

12. xy' - y = 2xln x => y' - (±)y = 21n x. 

v(x) = e~iv = e~ lnx = 1. (l)y' - (?) 2 y = =ln x 

s(i -y) = J lnx =► 5 -y = [lnx] 2 + C=>y = x[lnx] 2 +Cx 

e x g-x 

1 + e x y — (1 +e x ) • 

e~ x => (e x + l)y = e _x + C 

14. e -x dy + (e~ x y — 4x)dx = 0 => ^ + y = 4x e x => p(x) = 1 , v(x) = eJ ldx = e x =>■ e x ^ + y e x = 4x e 2x 
=> gj (y e x ) = 4x e 2x =$■ y e x = J 4x e 2x dx =>■ y e x = 2x e 2x — e 2x + C => y = 2x e x — e x + C e -x 

15. (x + 3y 2 ) dy + y dx = 0 => x dy + y dx = -3y 2 dy => ^ (xy) = -3y 2 dy => xy = -y 3 + C 

16. xdy + (3y — x _2 cos x) dx = 0 => y' + (?)y = x~ 3 cos x. Let v(y) = e = e 31nx = e lnx3 = x 3 . 

Then x 3 y ' + 3x 2 y = cos x and x 3 y = J cos x dx = sin x + C. So y = x~ 3 (sin x + C) 

17. ( x + l)gj + 2y = x => y' + (;^gy)y = ^y. Let v(x) = e^^+r dx = e 21n ( x+1) = e ln ( x+1 ) 2 = (x + l) 2 . 

So y'(x + l) 2 + (^tiy (x + l) 2 y = (^TT) (x + l) 2 =4> £ [y(x + l) 2 ] = x(x + 1) => y(x + l) 2 = fx(x + l)dx 
=> y(x +l) 2 = y + y+ C=^y=(x+ 1)~ 2 + y + cj . We have y(0) = 1 => 1 = C. So 

y = ( x + l) _2 (y + y + l) 

18. Xg|+2y = x 2 + l => y' + (?) y = x+ i.Letv(x) = ef(D dx = e lnx2 = x 2 . Sox 2 y' + 2xy = x 3 + x 

=* £( x2 y) = x 3 + x => x 2 y = g + y+ C=>y = g + g + y We have y(l) = 1 =*. 1 = 1 + C + ? => C = ±. 

e n v — x 2 | 1 , 1 _ x 4 + 2x 2 + 1 

^Oy-T + SI+2 - 4? 

19. ^ + 3x 2 y = x 2 . Let v(x) = e/ 3x ~ dx = e x \ So e x3 y' + 3x 2 e x3 y = x 2 e x3 => d d ^e x3 yj = x 2 e x3 => e x3 y = ?e x3 + C. 

We have y(0) = —1 =>■ e° 3 (— 1) = |e° 3 + C — l = j+ C=>C = — | and e x3 y = |e x3 — | => y = ? — |e~ x3 

20. xdy + (y - cos x)dx = 0 =>• xy' + y - cos x = 0 => y' + (?)y = Let v(x) = e^x dx = e lnx = x. 

So xy' + x(?)y = cos x => g-(xy) = cos x => xy = J cos xdx => xy = sin x + C. We have y(|) = 0 =>• (| )0 = 1 + C 
=> C = — 1 . So xy = — 1 + sin x => y = -1 g sinx 

21. xy' + (x — 2)y = 3x 3 e^ x y' + (^?)y = 3x 2 e~ x . Let v(x) = ei(^) dx = e x_21nx = g. So 

gy' + g (^)y = 3 => l (y • g) = 3 =*► y • g = 3x + C. We have y(l) = 0^0 = 3(1) + C =}► C = -3 
=>y-g=3x — 3=>y = xV x (3x - 3) 



13. (1 + e x )dy + (ye x + e X )dx = 0 => (1 + e x )y' + e x y = — e x => y' = 
V (x) = JttTST) = e We x +l) _ g x + y 

(e x + 1 )y ' + (e x + l)( I ggf)y = (fggy(e x + 1) => £[ (e x + l)y] = - 



e- x + C e- x + C 
y e x + 1 1 + e x 
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22. ydx + (3x-xy + 2)dy = 0=>| + 3 ^f^=0^| + |-x=-?^|+(^-l)x=-^ 
P(y) = ^ — 1 =>■ f P(y)dy = 31n y — y =>■ v(y) = e 31ny_y = y 3 e _y 

y 3 e~ y x' + y 3 e _y — 1 jx = — 2y 2 e^ y => y 3 e~ y x = J— 2y 2 e~ y dy = 2e~ y (y 2 + 2y + 2) + C 
=> y 3 = 2 ( y2 + 2y + 2 )+ Cey _ W e have y(2) = -1 => -1 = 2(i-2 + 2)+c e -' ^ c = _ 4e and 

3 ^ 2 (y 2 + 2y + 2) — 4e y+1 

y x 



23. To find the approximate values let y n = y n _i + (y n -i + cos x n _i)(0.1) with xo = 0, yo = 0, and 20 steps. Use a 
spreadsheet, graphing calculator, or CAS to obtain the values in the following table. 



X 


y 


X 


y 


0 


0 


1.1 


1.6241 


0.1 


0.1000 


1.2 


1.8319 


0.2 


0.2095 


1.3 


2.0513 


0.3 


0.3285 


1.4 


2.2832 


0.4 


0.4568 


1.5 


2.5285 


0.5 


0.5946 


1.6 


2.7884 


0.6 


0.7418 


1.7 


3.0643 


0.7 


0.8986 


1.8 


3.3579 


0.8 


1.0649 


1.9 


3.6709 


0.9 


1.2411 


2.0 


4.0057 


1.0 


1.4273 







24. To find the approximate values let y n = y n _i + (2 — y n „i)(2x n _i + 3) (0. 1 ) with xo = —3, yo = 1, and 20 steps. Use a 
spreadsheet, graphing calculator, or CAS to obtain the values in the following table. 



X 


y 


X 


y 


-3.0 


1.0000 


-1.9 


-5.3172 


-2.9 


0.7000 


-1.8 


-5.9026 


-2.8 


0.3360 


-1.7 


-6.3768 


-2.7 


-0.0966 


-1.6 


-6.7119 


-2.6 


-0.5998 


-1.5 


-6.8861 


-2.5 


-1.1718 


-1.4 


-6.8861 


-2.4 


-1.8062 


-1.3 


-6.7084 


-2.3 


-2.4913 


-1.2 


-6.3601 


-2.2 


-3.2099 


-1.1 


-5.8585 


-2.1 


-3.9393 


-1.0 


-5.2298 


-2.0 


-4.6520 







25. To estimate y(3), let y = y n _i + ( Xn x ‘ 1 j (0.05) with initial values xo = 0, yo = 1, and 60 steps. Use a spreadsheet, 
graphing calculator, or CAS to obtain y(3) « 0.8981. 

26. To estimate y(4), let z n = y n _i + ^ x ° -1 + (0.05) with initial values xo = 1, yo = 1, and 60 steps. Use a 

spreadsheet, graphing calculator, or CAS to obtain y(4) « 4.4974. 

27. Let y n = y n _i + ( e „ n _ 1+ y„ _ 1+ i ) (dx) with starting values xo = 0 and yo = 2, and steps of 0.1 and —0.1. Use a spreadsheet, 
programmable calculator, or CAS to generate the following graphs. 

(a) rr , 





i i i i 









[-0.2, 4.5] by [-2.5, 0.5] 
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(b) Note that we choose a small interval of x-values because the y-values decrease very rapidly and our calculator cannot 
handle the calculations for x < — 1. (This occurs because the analytic solution is y = —2 + ln(2 — e“ x ), which has an 
asymptote at x = —In 2 « 0.69. Obviously, the Euler approximations are misleading for x < —0.7.) 



I I 



[-1,0.2] by [-10, 2] 



28. Let y n = y n _i — (dx) with starting values xo = 0 and yo = 0, and steps of 0.1 and —0.1. Use a spreadsheet, 

programmable calculator, or CAS to generate the following graphs. 

(a) (b) 





29. x 1 1.2 1.4 1.6 1.8 2.0 

y -1 -0.8 -0.56 -0.28 004 (ET 

^=x=>dy = xdx=>y=^+C;x=l and y = - 1 
=>-l = i+ C=^C=-|=> y( exact) = y - \ 

=> y(2) = | i is the exact value. 




30. x 1 1.2 1.4 1.6 1.8 2.0 

y -1 -0.8 -0.6333 -0.4904 -0.3654 -0.2544 

g = i => dy = ^dx => y = ln|x| + C; x = 1 and y = -1 
=>— l=lnl+C=>C=— 1=> y(exact) = ln|x| — 1 
=> y(2) = In 2 — 1 « —0.3069 is the exact value. 




31. x 1 1.2 1.4 1.6 1.8 2.0 

y -1 -1.2 -0.488 -1.9046 -2.5141 -3.4192 

I=xy^f =xdx^ln|y| = f +C 
2 2 2 

=> y = eT +c = er • e c = Cie^ ; x = 1 and y = — 1 
=> — 1 = Cie 1/2 => Ci = — e 1/2 y(exact) = — e 1/2 • e? 

= — e ( x2-1 ) /2 y(2) = — e 3/2 w -4.4817 is the 

exact value. 
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33. ^ = y 2 - 1 => y' = (y + l)(y - 1). We have y' = 0 => (y + 1) = 0, (y — 1) = 0 => y = —1, 1. 

(a) Equilibrium points are — 1 (stable) and 1 (unstable) 

(b) y' = y 2 - 1 => y" = 2yy' => y" = 2y(y 2 - 1) = 2y(y + l)(y - 1). So y" = 0 =>■ y = 0, y = -1, y = 1. 



**,-1 v-l 




34. ^ = y - y 2 =t> y' = y(l - y). We have y' = 0 => y(l - y) = 0 => y = 0, 1 — y = 0 =>■ y = 0, 1. 

(a) The equilibrium points are 0 and 1 . So, 0 is unstable and 1 is stable. 

(b) Let — ► = increasing, < — = decreasing. 




y' = y - y 2 => y " = y' - 2yy' =>• y" = (y - y 2 ) - 2y(y - y 2 ) = y - y 2 - 2y 2 + 2y 3 => y" = 2y 3 - 3y 2 + y 
= y(2y 2 - 3y + 1) => y" = y(2y - l)(y - 1). So, y" = 0 => y = 0, 2y - 1 = 0, y - 1 = 0 =>• y = 0, y = 



y= 1. 

Let — ► = concave up, < — = concave down. 




35. (a) Force = Mass times Acceleration (Newton's Second Law) or F = ma. Let a = ^ ^ ^ = v^. Then 

ma = — mgR 2 s~ 2 =>■ a = — gR 2 s -2 => v^ = — gR 2 s -2 =>■ v dv = — gR 2 s _2 ds => fv dv = J— gR 2 s~ 2 ds 
=>• ^ + Ci =>• v 2 = ^ + 2Ci = ^ + C. When t = 0, v = v 0 and s = R => vg = + C 

=*► C = v 2 - 2gR => v 2 = + v 2 - 2gR 
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(b) If Vo = ^/2gR, then v 2 = =*> v = ^ since v > 0 if v 0 > i^^gR. Then ^ => ^/s ds = ^/2gR 2 dt 

=t> J s 1/2 ds = / i/2gR 2 dt => |s 3/2 = ^2gR 2 t + Ci =» s 3/2 = (§V2gR 2 )t + C; t = 0 and s = R 

=4> R 3/2 = (§^2gR 2 ) (0) + C => C = R 3/2 => s 3/2 = (fv/2gR 2 )t + R 3/2 = (§ R /2g)t + R 3/2 
= R 3 / 2 [(3R- 1 / 2 v ^) t+ i] = R 3 / 2 [(^p)t+l] = R 3 / 2 [(^)t+l] ^s = R [t + (^) t ] 2/3 

36. 'f' = coasting distance => (a86) [ 3a84) = 0.97 27.343. s(t) = ^(l - e - (k/m)t ) => s(t) = 0.97 (l - e - (27343/3a84)t ) 

=> s(t) = 0.97(1 — e -o.8866t^ a graph of the model is shown superimposed on a graph of the data. 
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t. (a) % = k£(c - y) => dy = -k£(y - c)dt ^ = -k^dt =>• f ^ = -f k^dt => ln|y - c| = -k£t + Ci 
=> y — c = ± e Cl e~ k v t . Apply the initial condition, y(0) = yo=^yo = c + C=>C = yo — c 
=> y = c + (y 0 - c)e _k v l . 

(b) Steady state solution: y^ = limy(t) = lim [ c + (y 0 — c)e~ k v l ] = c + (yo — c)(0) = c 



d(mv > = p + fv + ul — => F = d(mv) - (v + ul — =» F = m— 4- v— - v— - u— F = m— - u— 

dt r ^ v v ^ dt ^ r dt ! v ^ dt ^ r m dt T v dt v dt u dt ^ r m dt u dt ' 

^ = — b => m = — |b|t + C. At t = 0, m = mo, so C = mo and m = mo — |b|t. 



Thus, F = (mo - |b|t)^ - u|b| = -(m 0 - |b[t)|g| => £ = -g + 



=> V = -gt - ulnf + Ct 

mo-|bjt & \ m J 1 



v = 0att = 0=»Ci = 0. So v = — gt — u In 



t = 0 => y = -|gt 2 + c 



). Sov=-gt-uln(5^M) 

- t+ (^)in(^)- 



y=f — gt — uln^ m ° n J b ^ dt and u = c, y = 0 at 



3. (a) Let y be any function such that v(x)y = f v(x)Q(x) dx + C, v(x) = e / p ( x)dx . Then 

A(v(x) • y) = v(x) • y' + y • v'(x) = v(x)Q(x). We have v(x) = e/ p(x)dx => v'(x) = = e/ p{x)dx P(x) = v(x)P(x). 
Thus v(x) • y' + y • v(x) P(x) = v(x)Q(x) => y' + y P(x) = Q(x) => the given y is a solution. 

(b) If v and Q are continuous on [ a, b ] and x e (a, b), then d d [ J v(t)Q(t) dt = v(x)Q(x) 

=> J v(t)Q(t)dt = Jv(x)Q(x) dx. So C = yov(x 0 ) — Jv(x)Q(x) dx. From part (a), v(x)y = Jv(x)Q(x) dx + C. 
Substituting for C: v(x)y = J v(x)Q(x) dx + yov(xo) — J v(x)Q(x) dx =>• v(x)y = yov(xo) when x = xq. 



4. (a) y' + P(x)y = 0, y(xo) = 0. Use v(x) = e^ p ^ dx as an integrating factor. Then ^(v(x)y) = 0 =>■ v(x)y = C 

=> y = Ce~f p W dx and yi = Ci e~i p(x) dx , y 2 = C 2 e - i p ( x ) dx , yj(x 0 ) = y 2 (x 0 ) = 0, yj — y 2 = (Ci - C 2 )e“^ p 6) dx 
= C 3 e~/ p W dx and yi — y 2 = 0 — 0 = 0. So yi — y 2 is a solution to y ' + P(x)y = 0 with y(xo) = 0. 

(b) s(v(x)[yt(x)-y 2 (x)]) = ^(e/ p W dx [e-/ p W dx ( Cl -C 2 )]) = £(Ci - C 2 ) = £(C 3 ) = 0. 

f E( v ( x )[yi( x ) -yiWDdx = (v(x)[yi(x) — y 2 (x) ]) = f Odx = C 
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(c) yi = Ci e~/ p W d \ y 2 = C 2 e"/ p W dx , y = y t - y 2 . So y(x 0 ) = 0 => C 1 e~i’ p(x ) dx - C 2 e~f p W dx = 0 
=> Ci - C 2 = 0 => Ci = C 2 => yi(x) = y 2 (x) fora < x < b. 



5. (x 2 +y 2 )dx + xydy = 0 =>• g 



xy 



1 = — i 2 = pD) => F(v) = — i — v — 

x y/x x \x/ v / v x 



dx i dv 



dx 



dv 






ln|x 4 | + In 



2y 2 + x 2 



y/x x V x / ~ ^ V / v v x ^ v-F(v) 

0=>f f + / 2^1 =C=>ln|x| + iln|2v 2 +l|=C=^41n|x|+ln|2(^) 2 +l| =C 
= C => ln|x 2 (2y 2 + x 2 )| = C => x 2 (2y 2 + x 2 ) = e c =>• x 2 (2y 2 + x 2 ) = C 



= 0 



6- x=dy + (, 2 -x,)dx = 0^ S = ^4^g = -® 2 + l = FH)^F(v) = 
=*■/?+/* = c => ta M - ; = c => ln|x| - Jy = C => ln|x| 1 = C 



V — 



dv 

V- (— V 2 + v) 



= 0 



7. (xe y / x + y)dx - xdy = 0 => % = = e y / x + * = F( y ) =» F(v) = e v + v=>f + — 



dv 



e v + v) 



= o 



ff~f^= C => ln l x l + e“ v = C => ln|x| + e- y / x = C 



8. (x + y)dy+(x^y)dx = 0=> g = ^ = = F(|) ^ F(v) = =* f + - 



dv 






x + y 
dv 



= 0 



v 2 + 1 



+ J = 0 => ln|x| + tan 'v + iln|v 2 +1| = C 



21n|x|+2tan 'v + lnd) 2 -!-! = C => ln|x 2 | + 2tan 1 (|) + In y 1j x = C=>2tan '(D+ln 



y 2 +x 2 



= C 



dx j dv q 

x ' v — (v + cos(v — 1)) 



9. y' = l + cos(^) = \ + cos(f - 1) = F( y ) => F(v) = v + cos(v - 1) 

=> f y — f sec(v — 1) dv = 0 => In | x | — ln|sec(v — 1) + tan(v — 1)| = C => ln|x| — In I sec — l) + tan(| — l) 



10. (xsin| — ycosj) dx + xcos^ dy = 0 => ^ ^ S ‘xcos y = x — tan x = F(x) ^ F( v ) = v — tanv 

v — (/— tan v) = 0 => If + /cotvdv = 0 =* ln|x| + ln| sin v| = C =► ln|x| +ln| 



dx 



sin ; 



= c 
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